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A numerical model of a solar downdraft gasifier of biomass char (biochar) with steam 
based on the systems kinetics is developed. The model calculates the dynamic and steady 
state profiles, predicting the temperature and concentration profiles of gas and solid 
phases, based on the mass and heat balances. The Rosseland equation is used to calculate 
the radiative transfer within the bed. The char reactivity factor (CFR) is taken into account 
with an exponential variation. The bed heating dynamics as well as the steam velocity 
effects are tested. The model results are compared with different experimental results 
from a solar packed bed gasifier, and the temperature profile is compared to an experi¬ 
mental downdraft gasifier. Hydrogen is the principal product followed by carbon 
monoxide, the carbon dioxide production is small and the methane production is negli¬ 
gible, indicating a high quality syngas production. By applying the temperature gradient 
theory in the steam-only gasification process for a solar gasifier design, a solar downdraft 
gasifier improves the energy conversion efficiency by over 20% when compared to a solar 
packed bed gasifier. The model predictions are in good agreement with the experimental 
results found in the literature. 

© 2011 Elsevier Ltd. All rights reserved. 


1. Introduction 

The conventional autothermal gasification processes burn 
part of the carbonaceous compound in order to supply the 
energy necessary to enhance gasification reactions. Solar 
gasification seems to be a great solution in order to produce 
gaseous fuels, in which carbonaceous compounds are used 
exclusively as the carbon source, and the solar heat is used as 
the energy source for the endothermic reactions [1,2]. 

The fuel produced in this process is a high quality syngas, 
which is applicable for Fischer-Tropsch process or for power 
generation in fuel cells [3]. The most important advantage of 
solar biomass gasification is that they are available sources in 
wide areas, promoting energy independence and renewable 
energies [4], 


Using hydrogen as energy vector reduces carbon dioxide 
emissions and in a prospective way, the electricity can be 
stocked [5], Thus, the research has been concentrated in 
a hydrogen rich gas production. Some authors agree that 
when the gasification is done with steam the hydrogen yield 
increases (e.g. [6—11]). 

By the theoretical studies about gasification [12] it could be 
concluded that in order to improve the gasifiers performance, 
the cool phase should enter into the reactor at the side of the 
hottest point of the hot phase in order to enhance the endo¬ 
thermic reactions, while the exothermic reactions take place 
in the coolest points of the reactor, normally the output of the 
hot phase. 

This could be explained by the duality between the principal 
steam gasification reactions. The water-gas primary reaction is 
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Nomenclature 

A r Bed cross sectional area (m 2 ) 

A s Solid phase area in the cross sectional area (m 2 ) 

C ie Concentration ith gas in the emulsion phase 

(kg m- 3 ) 

C g Total gas concentration (kg mT 3 ) 

Cf Current value of concentration in the iteration 

(kg m- 3 ) 

Cp ep Specific heat of the emitter plate (kj kg -1 K _1 ) 

Cpi Specific heat of ith gas (kj kg -1 K _1 ) 

D Bed diameter (cm) 

D ie Diffusion of the ith in the emulsion phase (m 2 s _1 ) 
d p Particle diameter (m) 

g Acceleration of gravity (m s~ 2 ) 

AG? Free energy of formation of compound i (kj kg -1 ) 
H Bed height (m) 

h se Convection heat coefficient between solids and 

gas (kj s _1 nr 2 KT 1 ) 

AH r j Enthalpy of the jth reaction (kj kg -1 ) 

K Extinction coefficient (itt 1 ) 

Kg Thermal conductivity (W m~ a K _1 ) 

M Total solid mass in the reactor (kg) 

M ep Emitter plate mass (kg) 

m Mass flow rate (kg s _1 ) 

m e Biochar mass input flow rate (kg s _1 ) 


m s Biochar mass output flow rate (kg s _1 ) 

Q r Radiative flux density in the bed (W m~ 2 ) 

Or solar Concentrated thermal radiation (W nr 2 ) 

R Universal gas constant (kj mol -1 K _1 ) 

r ie Rate of the ith reaction in the emulsion gas 

(kg nr 3 s- 1 ) 

r s Rate of the ith reaction in the solids (kg m~ 3 s _1 ) 

t Time (s) 

T e Emulsion gas temperature (K) 

T ep Emitter plate temperature (K) 

T s Solids temperature (K) 

T 0 Input steam temperature (K) 

U Gas input superficial velocity (m s _1 ) 

Z Axial direction inside the reactor (m) 

Greek symbols 

a ij Stoichiometric coefficient of the component i in 

the j reaction (dimensionless) 
s Bed voidage 

e Solids emissivity 

Pg Gas density (kg nr 3 ) 

p s Solid density (kg nr 3 ) 

p Process efficiency 

a Stefan-Boltzmann constant (W nr 2 KT 4 ) 

Pg Gas viscosity (Pa s) 


endothermic, while the water shift reaction is exothermic. In 
the steam gasification process the first reaction takes place 
during the first contact between the steam and the solids; this 
is why it should be at the hottest point inside the reactor. 

On the other hand, when enough carbon monoxide is 
produced in order to react with the steam excess, the 
temperature of reactor should be lower in order to enhance 
the exothermic reaction, which is normally in the reactor 
output. This phenomenon could be named the temperature 
gradient theory in the steam-only gasification process. 

Different downdraft models have been proposed [4,13—19] 
as this kind of reactor has the advantage of satisfying, in 
general, the first contact condition whether the hottest phase 
is the gas or the solids. 

Giltrap et al. (2003) [14] introduced the concept of the char 
reactivity factor (CFR) which states the relative reactivity of 
the different chars; however this parameter failed to 
acknowledge the heat transfer in the reactor because it was 
taken as constant. Babu et al. (2006) [15] recognize the CFR as 
the key parameter in modelling a downdraft gasifier, and 
proposed to vary the CFR in different ways (constant, linear 
and exponential). They found that the linear and exponential 
variation gave the best predictions of the temperature and 
concentrations profiles compared to the experimental data 
reported in the literature. 

A solar packed bed gasifier has been proposed by Piat- 
kowsky et al. (2009) [3], in which the reactor consists of a 3D 
compound parabolic concentrator (CPC), two cavities sepa¬ 
rated by a SiC-coated graphite plate; with the upper one 
serving as the radiative absorber and the lower one containing 
the solids, the steam is injected at the reactor base. 


Dupont et al. (2007) [20] have done a time characteristic 
analysis of the steam gasification. They found that for 
temperatures between 1023 and 1273 K, the ratio between the 
chemical control and the external mass transfer control is in 
the order of 10 3 . This means that under these conditions the 
chemical regime would be the controlling step in the whole 
mass transfer. 

This paper aims to test the possible improvements of 
a solar packed bed gasifier performance by changing the 
reactor set-up to a downdraft reactor, in order to verify the 
temperature gradient theory in the steam-only gasification 
process. The model is based in the reaction kinetics with the 
CFR varying exponentially. The chemical regime is taken into 
account as the controlling phenomena for the mass transfer. 
The radiative heat transfer effects are included. 


2. Model development 

The solids are preheated to 473 K with an inert gas to dry them 
and to inhibit any eventual steam condensation. The model 
simulates the gasifying process of biochar. The pyrolysis and 
cracking reactions were not considered, as these two steps are 
supposed to take place in the preheating of the solids. 

The model uses the reactions kinetics proposed by Wang 
and Kinoshita (1993) [21], The reactions used for the model are 
shown in Table 1. 

The gasifier simulated is a downdraft reactor, in which an 
emitter plate heats the solids in the upper side of the reactor 
(see Fig. 1). The emitter plate is irradiated directly with 
concentrated solar energy. 
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Table 1 

— Considered chemical reactions. 


R n° 

Name 

Chemical reaction 

AH°298K 

(kj/mol) 

R1 

Water-gas 
shift reaction 

co + h 2 o«co 2 + h 2 

-41.1 

R2 

Water-gas 
(primary) reaction 

h 2 o + c«h 2 + co 

131.3 

R3 

Steam reforming 
reaction 

CH 4 + H 2 0«C0 + 3H 2 

206 

R4 

Boudouard 

reaction 

C0 2 + C«2C0 

172.8 


2.1. Reactor assumptions 

The following assumptions are made regarding the reactor 
operation: 

- No inert gas is used in the gasification process 

- Ideal behavior of gases is considered 

- The system parameters change only in the Z direction (see 

Fig- 2) 

- Char particles are spherical and of uniform size 

- The system is preheated with an inert gas to inhibit steam 
condensation, as well as to pyrolyse the biomass before the 
gasification process 

- Biochar contains only carbon 

- The heat transferred to the walls is not taken into account 


2.2. Reactions kinetics 

The system is completely represented by four stoichiometric 
independent reactions summarized in Table 2. 

When the gasification process is carried out in tempera¬ 
tures between 600 °C and 800 °C, the most important reactions 
are R1 and R2 because R3 needs high pressures and R4 needs 
high temperatures [12]. 

The char reactivity factor (CRF) is calculated exponentially 
as follows: 

CRF = e 00074Z 

Where Z is in mm. 

2.3. Heat and mass equations 

Fig. 2 shows a volume control (A R AZ ) fixed in the fluidized bed. 
The mass and heat balances are done for this volume control 
as follows: 

The generic mass balance of the ith gas including the 
chemical reactions will be: 

[Accumulation rate] = [Convective transfer] 

+ [Diffusion transfer] 

+ [Chemical generation] (1) 

Where the convective transfer is due to the velocity, the 
diffusion is described by the Fick law, the convection transfer 


Concentrated 
Thermal radiation 



Fig. 1 - A downdraft gasifier with concentrated thermal radiation as source of energy. 























BIOMASS AND BIOENERGY 35 (20Il) 2034-2043 


2037 


Z t 


AZ 


Z +AZ A 


u 


m 



Trans 

erGS 



<- 

-► 


Gas 

Solids 






-W- 


e (1—e) 

Fig. 2 — Volume control in the reactor. 


The generic heat balance of the ith gas, including the 
chemical reactions, are: 

[Accumulation rate] = [Heat flow input] - [Heat flow output] 

+ [Convection transfer between phases] 

+ [Chemical generation] (4) 


• Heat balance for the ith species in the gas phase 

Zl t (eCiC P i(T g )T g ) = C Pi (T g )T 3 j 

“ “ Ts ) + £ X/ r J AH rj ( 5 ) 

• Heat balance for the solids 

— £ )psCpsT s ) = Y f (^eCPsfs — WlsCPsT s) 

1 dA . 4 

+ A« ( 6 ) 

Where Qr is the radiative flux density, which is given by the 
Rosseland (1936) [22] approximation: 


between the phases is due to the concentration difference and 
the chemical generation is due to the chemical reactions. 


Qr(Z) = 


16 ffTg 9 T S 
3K 0Z 


( 7 ) 


• Mass equation for the ith species in the gas phase 


The initial and boundary conditions for the equation system 
are: 


8(eQ) _ 8 / 3CA 
8t 8Z\ 1 8Z ) 


V«& + e±c 


8Z 


j =i 


( 2 ) 


at t = 0 


T e = T 0 
T; = T S 0 


at Z = 0 and t > 0{T e = T 0 at Z = H and 


• Mass equation for the solids 


8M 

~dt 


dm 

8Z 



( 3 ) 


The initial and boundary conditions for the mass equations 
are: 


att = 0 


> 0 



at Z = 0 et t > 0{Qc = Q 0 at Z = H et t 


The condition at Z = H means that no further mass transfer 
is done towards the bottom of the reactor, then, it is an 
impervious surface. 



8T e 

8Z 

8^ 

8Z 


= 0 
= 0 


The condition at Z = H means that no further heat transfer 
is done towards the bottom of the reactor, then it is an isolated 
surface. 

For the temperature of the solid at Z = 0 (where the solids 
are irradiated) the following expression has been used for the 
boundary condition: 


^((1 - e)p s Cp s T s ) = -J- + (m s Cp s T s - m e Cp s T amh ) 

, 1 dAs, 1 dA ’ 

+ A„ dZ hse(Te Ts) A r dzg^ 


A r dZ 

+-( T t P - 7 t) 


( 8 ) 


Table 2 

- Kinetic parameters of reaction. 




R n" 

Rate of reaction 

Kinetic coefficients 

Units 

Reference 

R1 

H = CRF fci ^x co xh 2 o - X ^ XHi j 

k 4 = 2.824xl0- 2 e(- %P) 

kmol m~ 3 s' 1 

[22] 

R2 

Yi — CRF J?2 ^Xh 2 o ^ ^ 

k 2 = 1.517xl0 4 e( 1 ^ 2 ) 

kmol m' 3 s' 1 

[22] 

R3 

rv ( H 2 0 x \ 

r 3 = CRF k 3 (^XHjoXco - J 

k 3 = 7.310xl0- 2 e(fi5) 

kmol m' 3 s' 1 

[22] 

R4 

r 4 = CRF fe 4 

\ &4eq/ 

k 4 = 36.16e(-%P) 

kmol m' 3 s' 1 

[22] 
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The emitter plate temperature could be calculated as 
follows: 

df'ep _ 0/ solar /q\ 

dt MgpCpep 

The equilibrium constants in Table 2 are calculated as 
follows [17]: 




( 10 ) 


The NASA polynomials have been used to calculate the free 
energy of formation: 


AG° F' 

= A; + B,T + QT 2 + DjT 3 + EjT 4 + ^ + Giln(T) (11) 

Where A ; , B;, Q, D ; , E;, F; and G; are reported in [17]. 


3. Numerical solution 

The implicit volume finite method is used to estimate the 
solution of the equations system. The upwind method as 
described by Patankar (1980) [23] is also applied to the 
numerical solution. 


4. Results and discussion 

Two different gasifiers were simulated; the first simulates 
a solar downdraft gasifier with the dimensions reported by [3], 
the second simulates a solar downdraft gasifier with the 
optimum gasification length of the reduction zone reported by 
[16]. Table 3 shows the operational parameters simulated. The 
bed porosity and the bed emissivity are taken from [24]. 

Two emitter plate temperature profiles were tested order to 
study the influence of heating dynamics in the system. Fig. 3 
shows the emitter plate temperatures with time for the 
heating dynamics, these temperatures dynamics were taken 
from [3] for the cases of high carbon content feedstocks. The 
emitter plate temperature was used directly in equation (7). 

The heating dynamics 2 (HD2) heats the solids gradually 
until the emitter plate temperature reaches 1,700 K. On the 
other hand, the heating dynamics 1 (HD1) heats the solids 
faster to the same emitter plate temperature. For each heating 
dynamics, three different gas velocities were tested in G1 and 
they are listed in Table 3. 

The gas flow evolution with time for the three steam 
velocities for G1 and HD2 are shown in Fig. 4. There is not 
significant gas production before 20 min of gasification while 
the bed is heated. Between 20 and 40 min the gas production 
rises strongly as the solids temperature is increased (see 
Fig. 5). Then, as the solids temperature stabilizes the gases 
production slope gradually decreases and finally reaches the 
steady state. 

The principal gas produced is hydrogen followed by carbon 
monoxide, indicating a good syngas quality. Due to the fact 
that no combustion was conducted, the carbon dioxide yield is 
small for all runs. These trends are in good agreement with the 
results found experimentally by [3] for high carbon content 
feedstocks, where the hydrogen has the main concentration, 


Table 3 - Operational parameters. 


Parameter 

Gasifier 1 (Gl) 

Gasifier 2 (G2) 

Diameter (cm) 

15 

15 

Bed Height (cm) 

8.3 

23.6 

Steam input 

473 

473 

temperature (K) 

Steam velocities 

0.14-0.21-0.28 

0.69 

(m/s) 

Bed porosity 

0.5 

0.5 

Heat transfer 

0.054Re 1 ' 48 Kge/dp 

0.054Re 148 Kge/dp 

coefficient 

Bed emissivity 

0.75 

0.75 


the carbon monoxide concentration is bigger the carbon 
dioxide and the production of hydrocarbons is small. 

The gas flows for U = 0.1378 m/s are bigger than those 
reported by [3,12], indicating that a downdraft set-up could 
improve the gasifier performance compared to the packed and 
the fluidized beds, when the energy source is at the top of the 
reactor. This could be explained since the point of view of the 
temperature gradient theory, which states that a bigger 
temperature in the first contact point between the gas and the 
solids improve the hydrogen production. The fact of entering 
the gas at the top of the reactor, where the solids have the 
biggest temperature inside the reactor, ensure the best 
conditions of the first contact at the reactor input, thus 
increases the gas flows of the products. Fig. 6 shows the 
comparison of the present model results and those from 
Piatkowski et al. (2009) for high carbon content feedstocks. 

When the residence time of the steam is decreased (bigger 
steam velocities), the gas production decreases as well. There 
are two principal reasons for this: the time of the steam to 
react with the biochar is reduced and the heat transfer 
between the steam and the solids is improved (see Fig. 5), thus 
the temperature of the solids at the bed top is reduced and the 
yield of R2 is reduced as well. 

When the steam velocities are small, the bed heating is 
slower, this creates a bigger temperature gradient between the 
bed top and bottom before reaching the steady state, this 
gradient of temperature enhances the two principal steam 
gasification reactions in both extremes of the reactor, leading 
to a gas with higher hydrogen content at the reactor output, 
but at the same time a bigger carbon dioxide production. 



time (min) 


Fig. 3 - Emitter plate temperature profiles used in the 
simulations. 
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Fig. 4 - Gas flow evolution with time for Gl, HD2 (a) U = 0.1378 m/s, (b) U = 0.2067 m/s, (c) U = 0.2756 m/s. 
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Reactor height (m) time (min ) 



C 




Reactor height (m) 



Fig. 5 - Solid temperature in the bed evolution with the reactor height and time for Gl, HD2, (a) U = 0.1378 m/s, 
(b) U = 0.2067 m/s, (c) U = 0.2756 m/s. 
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South african coal Beech charcoal Present model 

Fig. 6 — Comparison of the present model higher flows and 
the experimental from Piatkowski et al. (2009) for high 
carbon content feedstocks (South African coal and Beech 
Charcoal). 


The energy conversion efficiency is calculated with the 
ratio between the energy content in the produced gas and the 
energy introduced to the system in steady state, as follows: 


V = 


Qsolar I filjeedstocfeLHVjeedstocfe 


( 12 ) 


The energy conversion efficiency values for the different 
runs are shown in Fig. 7. It is shown that the system efficiency 
is improved for a downdraft reactor, in which this parameter 
could be as high as 55% for small steam velocities compared to 
the packed bed where the efficiencies obtained by [3] for the 
high carbon content feedstocks are 23.3 and 29%. 

Fig. 8 shows the molar flow rates in steady state of 
hydrogen and carbon monoxide. The endothermic reactions 
could be closer to the equilibrium when the bed is heated 
gradually, and then bigger gas yields could be obtained, thus 
improving the process efficiency (see Fig. 7). 

Fig. 9 shows the solid temperature evolution with time. For 
small steam velocities the temperature gradient within the 
bed is significant from the start of the gasification. For the first 
20 min, while the temperature of the solids is under 800 °C at 
any point of the reactor, the gas production is low due to the 
insufficiency of energy to enhance R2. After 20 min the solid 
temperature at the bottom of the reactor begins to rise, at this 



■ U=0.1378 m/s, HD2 

■ U=0.2067 m/s, HD2 

■ U=0.2756 m/s, HD2 

■ U=0.1378 m/s, HD1 

■ U=0.2067 m/s, HD1 

■ U=0.2756 m/s, HD1 



Steam velocity (m/s) 


A CO HD1 
—■— H2 HD1 
■■■A- COHD2 
■ ■ ■■ ■ H2 HD2 


Fig. 8 - Molar flow rates of the principal gases in steady 
state. 


moment the gas production is at its maximum. After 50 min of 
gasification, the solid temperature gradient remains constant 
and close to 300 K until it reaches the steady state. These 
results are in good agreement with the experimental results 
found by [3]. 

A run was done in G2 in order to simulate the downdraft 
gasifier presented by Jayah et al. (2003) [16], Fig. 10 shows the 
temperature profiles from the present model and the experi¬ 
mental results obtained by Jayah et al. (2003) [16], The main 
goal of this comparison is not to validate the model results, 
but to check that, as reported by Babu et al. (2006) [15], an 
exponential variation of CFR is quite satisfactory and realistic 
to predict the temperature profiles, which is the normal 
behavior in a downdraft gasifier whether it is auto or alio 
thermal. A validation with these results is not relevant in this 
case, because Jayah et al presented experimental results for 
air-based autothermal gasification, which is not the same for 
steam-only endothermal gasification. 

A comparison of the gas yields is not possible because in 
the work of Jayah et al. (2003) [16] a combustion is conducted 
before the gasification zone. This changes the gases 



time (min) 


Fig. 7 - Energy conversion efficiency calculated with 
equation (11) for Gl. 


Fig. 9 - Temperature evolution with time at the reactor top 
and bottom for Gl, HD2 and U = 0.1378 m/s. 
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Fig. 10 - Temperatures profiles in steady state for G2 and U = 0.689 m/s. 


concentrations in the gasification zone input, thus the yields 
will be different compared to a steam-only gasification. 


5. Conclusions 

The development of a numerical model of a solar downdraft 
gasifier for gasifying biomass char (biochar) with steam based 
on the systems kinetics is presented. The model, based in the 
gasification kinetics, mass and energy balances, predicts gas 
yields and temperature profiles. The implicit volume finite 
method is used to estimate the solution of the equations 
system. 

The downdraft set-up could be a great solution in order to 
improve the performance of the packed bed and fluidized bed 
gasifiers with concentrated solar radiation in the upper side of 
the reactor. The gas produced is a high quality syngas, in 
which the hydrogen is the principal component followed by 
carbon monoxide; the carbon dioxide yield is small because no 
combustion is conducted. 

The system efficiency could be as high as 55% for small 
steam velocities. The energy conversion efficiency decreases 
when the steam velocity is increased and when the bed is 
heated quickly. 

The model predictions for the temperature profiles in G2 
are in very good agreement with the trends found experi¬ 
mentally and reported in the literature. Moreover, varying CRF 
exponentially improves the representation of the heat trans¬ 
fer throughout the bed. 

The influence of the walls in the heat transfer has not been 
taken into account, the distance between the the emitter plate 
and the top of the bed as well as at the bottom of the reactor is 
taken small enough to avoid this phenomena, this influence 
could be amplified if the distance is increased [25]. 

With the results reported in this paper, it is proved that taking 
into account the temperature gradient theory when designing 
the gasifier greatly improves the gasifiers performance. 
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